#define ddot_model(x) 3.918149e-08 *x*x*x+ -2.117145e-05 *x*x+ 2.318819e-03 *x+ 4.307391e-02
